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Abstract 

We study the Q 2 evolution of the nucleon's chiral-odd twist-3 distribution h^{x, Q 2 ) start- 
ing from the MIT bag model calculation. A simple GLAP equation for hi{x,Q 2 ) obtained 
at large N c is used for the Q 2 evolution. The correction due to the finite value of N c is 
0(1/N 2 ) ~ 10 % level. It turns out that the twist-3 contribution to hi(x, Q 2 ) is significantly 
reduced at Q 2 = 10 GeV 2 in contrast to the gi{x, Q 2 ) case. This is due to the fact that the 
corresponding anomalous dimension for hi, is larger than that for g 2 at small n (spin). 
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The EMC measurement of the nucleoli's g\ structure function |J inspired lots of theoret- 
ical activities on the nucleon's spin-structure functions in general as well as more precision 
measurements of g\ @. These structure functions provide us with a rich source of infor- 
mation about the spin distributions of quarks and gluons inside the nucleon. Jaffe and 
Ji || discussed general features of the quark distributions of the nucleon and relevant places 
where they can be measured. The nucleon has three independent twist-2 quark distributions, 
fi(x,Q 2 ) (spin-average), gi(x,Q 2 ) (helicity asymmetry), hi(x,Q 2 ) (helicity flip), and three 
independent twist-3 quark distributions e(x,Q 2 ), g 2 (x,Q 2 ), hi^x.Q 2 ). Twist-2 distributions 
have a simple parton model interpretation and contribute to various hard processes in the 
leading order with respect to (Q is the hard momentum of the external hard probe.) 

On the other hand, the twist-3 distributions represent complicated quark-gluon correlations 
in the nucleon, and is generally difficult to be measured, since they are often hidden behind 
the leading twist-2 contributions. However, g 2 and can be measured in the absence of the 
leading twist-2 contributions through the proper asymmetries in the polarized deep inelastic 
scattering and the polarized Drell-Yan process, respectively @, [|. In this sence, they are in- 
teresting higher twist distribution functions. In fact, E143 collaboration || presented a first 
nonzero data for g 2 , which anticipates a forthcoming significant progress in twist-3 physics. 

So far accumulated experimental data on f\ and gi allowed us to parametrize in the 
next-to-leading order for fi |J and for gi 0. But nothing is known about the actual shape 
of h%, g 2 and hi except some guess by the bag model calculations^, ^|, (Since there 
is no practical way of isolating e, it will not be considered in this work.) The bag model 



has been reasonably successful in describing various properties of hadrons ||10|| , and has been 



applied to calculate the structure functions ||, [9|, [11], [T^, [L3], |14|]. Since the bag model is a low 
energy effective hadron model, its prediction for the structure functions has to be evolved to 
higher scale to confront experimental data. After the Q 2 evolution, it could approximately 
reproduce the valence parts of fi and g\. The purpose of this short note is to present the first 
and a rough estimate of the magnitude of Hl{x, Q 2 ) at high Q 2 staring from the bag model 
calculation. We are especially interested in the speed of the Q 2 -evolution of h L compared 
with that of g 2 and the chiral-odd twist-2 distribution h±. Since it is not our purpose here 
to construct a more realistic model, we shall not persue the projection method to restore 
the translational invariance as was tried in |9|, 12 , [14| . We refer those attempts to future 



studies. 

We first recall the definition of the quark distributions in our interest 0] : 

^e^(PS\mi»l^{\n)\ Q \PS) 
In 

= 2 [ 9l (x, Q 2 ) PlM (S ■ n) + g T (x, Q 2 )S ± ^ + M 2 g 3 (x, Q 2 )S ■ nn^ 
^e tXx (PS\ij(0)a, u i l5 ^(\n)\ Q \PS) = 2 Ih^x, Q 2 ) (S XfiPu - S XuPfi ) /M 

+h L (x, Q 2 )M ( Ptl n u - p v n,j) (S ■ n) + h 3 (x, Q 2 )M {S^n v - S^n^) , (2) 
where \PS) denotes the nucleon (mass M) state with the four momentum P and the spin 
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S, and the two light-like vectors, p and n, are introduced by the relation P M = + ^n 11 , 
p ■ n = 1, p 2 = n 2 = 0. For the nucleon moving in the ^-direction, p = -^(1,0,0, 1) and 

'1,0,0,— 1). V — > oo corresponds to the infinite momentum frame and V = M/\/2 



n 



V2V 
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corresponds to the nucleoli's rest frame. is decomposed as = (S '-ri)jf + (S ■p)n 11 + S±. In 
([!]) and (H), lightcone gauge, n-A ~ A + = 0, was employed. The above distribution functions 
9i,t {qt — 9\ + g 2 ) and /i^x, etc are defined for each quark flavor ip = ip a (a = u, d, s, ...) and 
have support — 1 < x < 1 [|L5[]. The replacement ?/> a — > CY> a , ?/> a — > — ^ a C~ x defines the 
anti-quark distributions gf T (x) etc for each quark distribution g± T {x) etc. They are related 
as 9it(~ x ) = 9it{ x )i x ) = ~ F° r the polarized deep inelastic scattering, 

physically measurable structure functions are the combination X) Q ^(^i r( x ) +5 i it( x )) with 
the Bjorken x (0 < x < 1) and the electric charge of a (anti-)quark flavor a, e a . Here and 
below, we often suppress the explicit Q 2 dependence of the distributions. 

The Q 2 dependence of these structure functions is calculable in perturbative QCD. The 
twist-2 distributions, g\ and hi, obey simple GLAP equations |T6|. On the other hand, 
the Q 2 dependence of the twist-3 distributions, g 2 and h L , is quite sophisticated because 
the number of quark-gluon-quark operators increases with the moments (or spin). The 
calculation of the one-loop anomalous dimension matrix for all the twist-3 distributions has 
been completed |17|, [18|, [19|, ^(J, ^TJ, and an analogue of GLAP equation relevant to describe 
Q 2 evolution of the whole x dependent distributions has also been derived in ]TB|, |21[ . These 
equations for the twist-3 distributions are the evolution equation for the corresponding parent 
distributions and is not convenient for practical applications. However, there is a very useful 
news for physicists working on higher twist effects. It has been proved that at large -/V c , the Q 2 
evolution of all the twist-3 distributions is described by simple GLAP equations with slightly 
different forms for the anomalous dimensions from the twist-2 distributions p2|, |23"| , pCf| : [] 
The Q 2 evolution (for g 2 , only for nonsinglet piece) is given by 



M n [HQ 2 

M n \h L (Q 2 ) 



M, 



e(Q 2 



L^M n \~h L ^ 2 ) 



(3) 
(4) 
(5) 



where M n [g(Q 2 )] = j\dx x n g(x,Q 2 ), L = a s (Q 2 ) /a s (n 2 ), b = ^N c - £Nf. 



fN c - lN f . g 2 and h L 
denote the twist-3 parts of g 2 and hL, respectively. The corresponding anomalous dimensions 
in (H)- (H) are given by 




7 : = 2N c \S n - 7 



4 2(n+l) 



(6) 
(7) 
(8) 



with S n = E ■ 



Furthermore, these anomalous dimensions are the lowest eigenvalues of 



the anomalous dimension matrices at large N c . Since these relations are obtained by a mere 
replacement Cf = {N 2 — 1)/2N C — > N c /2 in the complete one-loop anomalous dimension 

1 In a recent work|24|], it has also been shown that the same simplification at large iV c occurs for the Q 2 
dependence of all the twist-3 fragmentation functions. 
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matrices at finite N c , the correction due to the finite value of N c is 0(1/N 2 ) ~ 10% level, 
which is sufficient for practical application. The essential ingredient in (ff)-(|8D is that a 
knowlegde on g 2 {x), h^ix) and e(x) at one scale is sufficient to predict them at an arbitrary 
scale, which is not the case at finite N c . This fact provides us with a useful framework to 
confront experimental data at various Q 2 of the twist-3 distribution. In fact (^) and (||) were 
used to predict the shape of g 2 at high Q 2 starting from the bag model calculation 0. A more 
favorable feature of hi compared with g 2 is that hi does not mix with gluon distributions 
owing to its chiral-odd nature. Therefore Q 2 evolution for h^ and e is given by (|[), (|[), (|7]) 
and (§) even for the flavor singlet piece and thus we can get more reliable and accurate form 
in the small x region compared to g 2 . This work is devoted to the study of Q 2 evolution of 
h L with (D and <g). 

In the rest frame of the nucleon, one can conveniently calculate the above distributions 
using the MIT bag model. The result for hi and hi with one quark flavor in the nucleon is 
given as M 



h L (x) 



ujMR 



2ti(uj - J\y min \ 



Here ti is given by 



ydy 



ujMR 

2ti(uJ - l)joH J\Vmin\ 



ydy 



t (uj,y) 2 - ti(u,yY 



JJr, 



+ t 1 (u,y) 2 [ V -™ 



(9) 
(10) 



ti(u,y) 



dz z 2 ji(uz)ji(yz), 



'ID 



where ji is the l-th order spherical Bessel function, and to is determined by the relation 
tana; = —uj/{uj — 1). For the lowest energy mode, uj = 2.04. y min is defined as y m in = 
MRx — uj with the bag radius R determined by the relation MR = Au. is decomposed 
into the twist-2 piece which can be expressed in terms of hi and a purely twist-3 piece hi as 



h L (x) 



2x / dy 



h(y) 



y l 



+ h L (x) 



-2x 



-i 



dy — 5— + h L {x). 

y 



< x < 1 
-1 < x < 



'121 



The bag model prediction above has to be regarded as a distribution at some low energy 



scale Q 2 = nl ag < 1 GeV 2 . For hi, we regard 
scale. 



as a valence distribution at this low energy 



In order to evolve the bag model prediction for hL from [i 2 ag to Q 2 according to (|]), we 



used a method in ||25|| . For the moment, we symbolically represent hi^(x) by h(x). If one 
defines h±(x) = h(x)±h(—x) = h(x)^fh(x), the even (odd) moments of h + (x) (h_(x)) on the 
interval [0, 1] agree with «M n [/i], whose Q 2 evolution is given by (f|) and (0) and its analogue. 
We assume Q 2 evolution of all the moments of h + (x) and h_(x) on [0, 1] is described by 
the same anomalous dimensions as was often assumed to describe Q 2 evolution of fi(x, Q 2 ) 
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and gi (x,Q 2 ), and construct h(x,Q 2 ) on [—1,1]. This is equivalent to assume that the Q 2 
dependence of the moments of h(x,Q 2 ) on [0, 1] and [—1,0] are separately governed by the 
same anomalous dimension in (^), which is a sufficient condition to satisfy (01). To describe 
the method in | 25fl , we introduce Bernstein polynomial defined by 



k\ 



i\(N-k-iy: 



and note that it satisfies the relation 



lim b { - N ' k \y) = 5{x-y) 

N,k— >oo 
k/N^x 

for < x,y < 1. Using (|4]), (|) and (0), we get 



(13) 



(14) 



h L (x,Q z 



,. (iV + i)!^ fc 
hm — > , 



N,k^oo 
k/N^x 



k\ 



1=0 



l\(N-k-l)\ Jo 



dyy k+l h L (y,Q 2 



(jV + l)!^ fc 
hm — > , 



N,k^oo 
k/N^x 



-Z^V 6 ° A 1 dyy k+l h L (y,n 2 ). (15) 



i=0 



l\(N — k — l)\ 



Since the summation over / in ([15]) oscillates due to the factor (—1)', the direct use of (|T5| ) 
is not convenient. To avoid this difficulty we shall utilize the following procedure. Expand 



L^/"° as 



Q(L) 



() (n + p) 4 +p( L )' 

where a(L), Ci(L) and p are the constants determined below. Then (|1^) is rewritten as 



(16) 



h(x,Q 2 ) 



—b{x,y;Q,ii)h{y,v), 

y 



b{x,y;Q,fi) = a(L) ( - ] X](ln- 



2/, 



i=0 



x 



where we have used the relation 
(jV + l)!^ fc (-1)' 



,k+l 



r(* + p) 



fcfl M p -V ln .y + -\ (18) 



(17) 



lim 

N,k~^oo 
k/N^x 



k\ 



§ V.(N-k-l)\ {k + l + py+p 



The expansion ([TB]) can be obtained by applying the following asymptotic expansion to 7^ 
in©: 



n+l 



7s + ln(n + 1) + 



2(n + 



T>-£ 



2/,- 



\ 2k(n + l) 2k 



(19) 
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where 7^ = 0.577... is the Euler constant and life's are the Bernoulli numbers. This proce- 
dure gives p = 1 and the coefficients in (fllf) . (See [^, for the details.) We have used 



first four terms (i.e. up to B s ) in the expansion (|T9|), which gives enough precision. 

Next we need to determine the bag scale fi bag . Phenomenological values for n 2 ag adopted 

,2 



in the previous studies scatters below 1 GeV : Jaffe and Ross [|26| extracted fi bag = 0.75 



GeV 2 from the 5-th moment of a F 3 data. Schreiber et al. Jl3| took fi? = 0.25 GeV 2 to 



reproduce approximate shape of F2 at Q 2 = 10 GeV 2 . More recently, Stratmann determined 
f^bag = 0.081 GeV 2 by comparing the second moment (momentum sum rule) of the bag 
model prediction with that of the valence distribution determined by Gliick et al. 0. f] Our 
purpose here is to see how evolves compared with other distributions such as hi and g2, 
and therefore we shall show the results for two values of jJ^ ag , ^ 2 ag = 0.081,0.25 GeV 2 , for 
future references. For other parameters, we set Nf = 3, Aq CD = 0.232 GeV in a s (Q 2 ). 

Figure 1 (a) and (b) show the results for xhi(x,Q 2 ) and xgi(x, Q 2 ) at Q 2 = 10 GeV 2 
with two values of /i 2 together with the bag calculations. For the Q 2 evolution, we have 



used the anomalous dimension for hi calculated in |27|, gi(x,Q 2 ) at Q 2 = 10 GeV 



is 



strongly peaked in the small x region, since the anomalous dimension for -Mol^i] is zero 
while it is 4/3 for A^o[^i]- (If one plots hi and gi instead of xhi and xg%, this feature is 
more conspicuous.) Figure 2 shows the bag calculation for hi decomposed into the twist-2 
and -3 contributions § ((a)) and h L evolved to Q 2 = 10 GeV 2 with Lt 2 bag = 0.081 GeV 2 ((b)) 
and fil ag = 0.25 GeV 2 ((c)). In Fig. 3, we plot only the twist-3 piece of hi, hi, taken 
from Fig. 2 to see how it evolves with Q 2 . These graphs show clearly that at higher Q 2 the 
contribution from hi(x, Q 2 ) is significantly reduced and h^x, Q 2 ) is dominated by the twist- 
2 contribution. Although our calculation starts from the bag model prediction, this tendency 
can be taken as model independent. Comparison of Fig. 3 and Fig. 1 shows that h^ evolves 



faster than hi as is expected from the magnitudes of the anomalous dimensions [|T9[ . For 
comparison we have also shown g2(x,Q 2 ) in Fig. 4. Since this distribution is accessible in 
the polarized DIS, we plotted the combination g2(x, Q 2 ) +g2{x, Q 2 ) = g2(x, Q 2 ) +g 2 (—x, Q 2 ) 
for one quark-flavor with unit charge. At the bag scale, twist-3 contribution g 2 is comparable 
to the twist-2 contributin as in the case of h^. This feature more or less survives even at 
Q 2 = 10 GeV 2 in contrast to the h L case. This is because 7^ > 7^ especially at small n and 
hence Q 2 evolution of hL in the small x region is faster than that of g2- As was stated before, 
flavor- singlet part of g2 mixes with the gluon distribution and Q 2 evolution of singlet g2 is 
not given by (|3|). Singlet ^2 is probably more enhanced at small x region. If the bag model 
gives a good description even for the twist-3 distribution h L , our present study indicates 
that it will be extremely difficulte to extract hL{x,Q 2 ) at high Q 2 . On the other hand, if 
future experiments show hi{x,Q 2 ) is still sizable at high Q 2 , it means that the naive bag 
model calculation is not suitable to describe quark-gluon correlation represented by hL in 
the nucleon. In any case, it is very interesting to confirm these general features in the future 
collider experiments. 



Comparison with other parton distributions in Q gives almost the same numbers for 



^bag • 



6 



References 



[1] EMC collaboration, J. Ashman et ai, Phys. Lett. B206 (1988) 364. 

[2] SMC collaboration, D. Adams et ai, Phys. Lett. B329 (1994) 399; B357 (1995) 248; 
E143 collaboration, K. Abe et a/., Phys. Rev. Lett. 74 (1995) 346; 75 (1995) 25; Phys. 
Lett. B364 (1995) 61. 

[3] R.L. Jaffe and X. Ji, Nucl. Phys. B375 (1992) 527. 

[4] R.L. Jaffe, Comm. Nucl. Part. Phys. 19 (1990) 239. 

[5] E143 collaboration, K. Abe et a/., Phys. Rev. Lett. 76 (1996) 587. 

[6] M. Gliick, E. Reya and A. Vogt, Z. Phys. C53 (1992) 127; A.D. Martin, W.J. Stirling 
and R.G. Roberts, Phys. Lett. B354 (1995) 155; H.L. Lai, et. al. Phys. Rev. D51 (1995) 
4763. 

[7] M. Gliick, E. Reya, M. Stratmann and W. Vogelsang, Phys. Rev. D53 (1996) 4775; 
T. Gehrmann and W.J. Stirling, Phys. Rev. D53 (1996) 6100; G. Altarelli, R. Ball, S. 
Forte and G. Ridolfi, |hep-ph/ 9701289] . 

[8] R.L. Jaffe and X. Ji, Phys. Rev. D43 (1991) 724. 

[9] M. Stratmann, Z. Phys. C60 (1993) 763. 
[10] A. Chodos et ai, Phys. Rev. D9 (1974) 3471. 
[11] R.L. Jaffe, Phys. Rev. Dll (1975) 1953. 

[12] C.J. Benesh and G. Miller, Phys. Rev. D36 (1987) 1344; D38 (1988) 48. 
[13] A.W. Schreiber, AW. Thomas and J.J. Londergan, Phys. Rev. D42 (1990) 2226. 
[14] A.W. Schreiber, A.I. Signal and A.W. Thomas, Phys. Rev. D44 (1991) 2653. 
[15] R.L. Jaffe, Nucl. Phys. B229 (1983) 205. 

[16] V.N. Gribov and L.N. Lipatov, Sov. J. Nucl. Phys. 15 (1972) 438; 
L.N. Lipatov, ibid. 20 (1975) 94; 

G. Altarelli and G. Parisi, Nucl. Phys. B126 (1977) 298. 

[17] E.V. Shuryak and A.I. Vainshtein, Nucl. Phys. B201 (1982) 141; 
PC Ratcliffe, Nucl. Phys. B264 (1986) 493; 
LI. Balitsky and V.M. Braun, Nucl. Phys. B311 (1988/89) 541; 
X. Ji and C. Chou, Phys. Rev D42 (1990) 3637; 
J. Kodaira, Y. Yasui and T. Uematsu, Phys. Lett. B344 (1995) 348; 
J. Kodaira, Y. Yasui, K. Tanaka and T. Uematsu, Phys. Lett. B387 (1996) 855. 

[18] A.P Bukhvostov, E.A. Kuraev and L.N. Lipatov, Sov. Phys. JETP 60 (1984) 22. 

[19] Y. Koike and K. Tanaka, Phys. Rev. D51 (1995) 6125. 

[20] Y. Koike and N. Nishiyama, Phys. Rev. D55 (1997) 3068. 

[21] A.V. Belitsky and D. Miiller, |hep-ph/9702354] . 

[22] A. Ali, V.M. Braun and G. Hiller, Phys. Lett. B266 (1991) 117. 

[23] LI. Balitsky, V.M. Braun, Y. Koike and K. Tanaka, Phys. Rev. Lett. 77 (1996) 3078. 



7 



[24] A.V. Belitsky, |hep-ph/9702356[ 

[25] A. Gonzales-Arroyo, C. Lopez and F. J. Yndurain, Nucl. Phys. B153 (1979) 161. 

[26] R.L. Jaffe and G.G. Ross, Phys. Lett. B93 (1980) 313. 

[27] X. Artru and M. Mekhfi, Z. Phys. C45 (1990) 669. 



8 



Figure captions 



Fig. 1 (a) Bag model calculation for xh±(x) (dash-dot line) and hi(x,Q 2 = 10 GeV 2 ) with 
nl ag =0.08l GeV 2 (solid line) and 0.25 GeV 2 (dashed line), (b) The same as (a) but for for 
xg^x.Q 2 ). 

Fig. 2 (a) Bag model calculation for h L (x) (solid line) decomposed into the twist-2 (dashed 
line) and the twist-3 (dash-dot line) contributions, (b) hi{x,Q 2 ) at Q 2 = 10 GeV 2 with 
f^lag = 0.081 GeV 2 decomposed into the twist-2 and twist-3 contributions, (c) The same as 
(b) but with n 2 bag = 0.25 GeV 2 . 

Fig. 3 Twist-3 contribution to Iil(x, Q 2 ) at the bag scale (dash-dot line) and Q 2 = 10 GeV 2 
with fx 2 bag =0.081 GeV 2 (solid line) and 0.25 GeV 2 (dashed line). 

Fig. 4 (a) Bag model calculation for g 2 (x) decomposed into the twist-2 (g^ W ( x )) an d the 
twist-3 (g 2 (x)) contributions, (b) g 2 (x,Q 2 ) at Q 2 = 10 GeV 2 with fi 2 bag = 0.081 GeV 2 . (c) 
The same as (b) but with (j% = 0.25 GeV 2 . 



9 




X 

Fig.1(a) 




Fig.3 



15 



